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Abstract.  Hartree-Fock  and  Kohn-Sham  orbital  energies,  the  latter  computed  with 
several  different  exchange/correlation  functionals,  are  compared  and  analyzed  for  twelve 
molecules.  The  Kohn-Sham  energies  differ  significantly  from  experimental  ionization 
energies,  but  by  amounts  that  are,  for  a  given  molecule  and  exchange/correlation 
functional,  roughly  the  same  for  all  of  the  valence  orbitals.  With  the  exchange/correlation 
functionals  used,  the  energy  of  the  highest  occupied  Kohn-Sham  orbital  does  not 
approximate  the  corresponding  ionization  potential  any  better  than  do  the  other  orbital 
energies. 

Key  words:  Hartree-Fock  orbital  energies,  Kohn-Sham  orbital  energies,  ionization 
potentials. 


Introduction 


The  Hartree-Fock  method  and  the  Kohn-Sham  version  of  density  functional  theory  have  in 
common  the  feature  that  they  both  involve  the  iterative  solution  of  a  set  of  one-particle 
eigenequations  [1-3],  The  resulting  eigenfunctions  { (|)i}  and  their  associated  eigenvalues 
{ Ei}  are  generally  designated  as  atomic  or  molecular  orbitals  (depending  upon  the  system 
in  question)  and  their  corresponding  orbital  energies. 

In  Hartree-Fock  theory,  the  physical  interpretation  of  these  orbital  energies  is  based 
upon  the  approximation, 

Ij  =  Ejjp(N-l,i)  -  Ejjp(N)  (1) 

in  which  li  is  the  ionization  potential  of  an  electron  in  orbital  (l)i,  Ehf(N)  is  the  energy  of 
the  N-electron  parent  system  and  EHF(N-l,i)  is  the  energy  after  removal  of  an  electron 
from  <|)i,  leaving  the  geometry  and  other  orbitals  unchanged.  Since  it  can  easily  be  shown 

that  [4,  5], 


EHF(N-l,i)  -  Ehf(N)  =  -eHF.i 

then  it  can  be  inferred  that  -Ehf.I  =  magnitudes  of  the  Hartree-Fock  orbital 

energies,  {lEnF,!'}  >  ®2n  be  viewed  as  approximately  equal  to  the  ionization  potentials,  (li), 

of  the  respective  electrons  [4, 5].  Some  justification  for  eq.  (1)  comes  from  Koopmans’ 
theorem  [4,  6],  which  assures  that  the  wave  function  obtained  by  simply  removing  an 
electron  from  (t)i  is  stable  with  respect  to  further  variation  in  (j)!.  However  eq.  (1)  still 
neglects  the  fact  that  loss  of  an  electron  produces  a  spatial  rearrangement  of  the  charge  in 
the  remaining  orbitals,  stabilizing  the  ion;  this  factor  in  itself  would  cause  I  eHF,il  to 
overestimate  li.  On  the  other  hand,  electronic  correlation,  which  is  omitted  in  Hartree-Fock 


theory,  is  expected  to  lower  Ehf(N)  more  than  EHF(N-l,i);  this  alone  would  make  lepjp^il 
underestimate  li.  In  Hartree-Fock  calculations,  therefore,  the  effects  of  these  two  sources 
of  error  should  peirtially  cancel  [2, 5]. 

There  has  been  considerable  discussion  of  the  physical  interpretation  (if  any)  to  be 
given  the  orbital  energies  { eKS.i}  in  Kohn-Sham  theory  [2, 3, 7-11],  in  which  Koopmans’ 
theorem  is  not  applicable.  The  general  conclusion  has  been  that  the  energy  of  the  highest 
occupied  orbital  (hoo)  in  exact  Kohn-Sham  theory  should  be  equal  in  magnitude  to  the  first 
ionization  potential, 

“^KS,hoo  ~  i^min 

but  that  no  specific  physical  meaning  can  be  assigned  to  the  other  orbital  energies  [2, 3, 7- 
12].  Eq.  (3)  requires  that  the  Kohn-Sham  effective  potential  approach  zero  asymptotically, 
a  point  that  has  recently  been  discussed  [12-14]. 

Our  present  objective  is  to  provide  and  to  analyze  quantitative  data  relevant  to  these 
points.  Our  primary  focus  will  be  on  approximate  Kohn-Sham  orbital  energies,  obtained 
using  several  different  exchange/correlation  functionals,  but  Hartree-Fock  values  are  also 
included,  for  completeness  and  for  perspective.  These  computed  results  will  be  compared 
to  experimentally-determined  electronic  ionization  potentials  for  a  group  of  twelve 
molecules. 

Procedures 

Hartree-Fock  and  Kohn-Sham  geometry  optimizations  were  carried  out  with  Gaussian  94 
[15]  for  eleven  polyatomic  molecules.  Four  different  basis  sets  were  used  for  the  Hartree- 
Fock  calculations:  STO-3G,  3-21G,  6-31G*  and  6-31-t-G**.  The  Kohn-Sham  results 
were  obtained  with  several  combinations  of  exchange  and  correlation  functionals,  including 
the  Becke  (B)  [16],  the  Becke  three-parameter  hybrid  (B3)  [17],  the  Perdew  86  (P86)  [18], 
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the  Lee-Yang-Parr  (LYP)  [19]  and  the  Perdew-Wang  91  (PW91)  [20].  The  basis  set  for 
most  of  the  Kohn-Sham  work  was  the  6-3 1+G**;  however  the  B3P86  calculations  were 
also  carried  out  at  the  6-31 1++G(3df,  3pd)  level,  to  assess  the  effect  of  a  rather  large  basis 
set.  The  computed  orbital  energies  are  compared  to  experimentally-determined  vertical 
ionization  potentials  obtained  through  photoelectron  spectroscopy;  the  latter  data  are  taken 
primarily  from  the  compilation  by  Duffy  and  Chong  [21]. 

Results 

Table  1  summarizes  our  results  for  the  eleven  molecules  treated  in  the  present  study  plus 
one  (N2O)  for  which  the  calculations  were  carried  out  earlier  [22].  For  each  molecule  are 
given  the  differences  between  the  experimental  ionization  potentials  [21]  and  the 
magnitudes  of  the  orbital  energies  as  computed  by  nine  different  Hartree-Fock  and  Kohn- 
Sham  procedures.  The  highest  occupied  orbitals  are  listed  first,  followed  by  several 
valence  orbitals  in  order  of  increasing  If.  The  calculated  total  molecular  energies  decrease 
froin  left  to  right  across  the  table,  the  one  exception  being  that  E(B3PW91)  <  E(BLYP)  for 
cyclopropene. 

Discussion 

For  the  twelve  molecules  in  Table  1,  the  magnitudes  of  the  Hartree-Fock  orbital 
energies  are  overall  closer  to  the  corresponding  ionization  potentials  than  are  the 
approximate  Kohn-Sham  energies.  The  Hartree-Fock  (except  for  the  STO-3G)  usually 
overestimate  and  the  Kohn-Sham  underestimate  the  experimental  values.  Thus  in  general. 


-eKS,i  <  li  <  -eHF,i 


(4) 


as  was  found  earlier  for  atoms  [8, 1 1]  and  for  N2O  [22]. 

The  Hartree-Fock  6-3 IG*  and  6-3 1+G**  results  are  quite  similar,  as  are  also,  for 
the  most  part,  the  3-2 IG.  These  three  basis  sets  are  all  relatively  effective  in  predicting  the 
ionization  potential  for  the  highest  occupied  orbital;  the  average  absolute  error  is  0.7  ev, 
and  would  be  0.5  ev  if  not  for  the  OF2  molecule.  (The  3-21G  basis  set  usually  gives  the 
best  agreement.)  The  predictions  tend  to  be  considerably  poorer  for  the  lower  orbitals,  but 
no  consistent  pattern  is  apparent.  (The  STC)-3G  results  are  particularly  erratic.)  For  a 
given  molecule  and  basis  set,  the  quantities  (-enF.i  “  li)  cover  rather  large  ranges, 

typically  2  -  3  ev. 

In  contrast  to  the  Hartree-Fock,  the  Kohn-Sham  results  do  show  certain  patterns: 

(1)  The  Kohn-Sham  orbital  energies  differ  from  the  ionization  potentials  by  at  least  2  ev 
and  often  much  more.  This  apparently  reflects  a  need  to  more  effectively  remove 
electronic  self-interaction  [8, 11].  However  the  quantities  (-£^51  -  Ij) ,  for  a  given 

molecule  and  exchange/correlation  combination,  show  a  notable  uniformity,  usually 
being  within  a  range  of  0.3  -  0.7  ev.  At  these  levels  of  approximation,  the  energy  of 
the  highest  occupied  orbital  does  not  provide  a  better  estimate  of  the  corresponding 
ionization  potential  than  do  the  other  orbital  energies.  This  is  true  even  when  a  large 
basis  set,  6-311++G(3df,3pd),  is  used. 

(2)  The  closest  agreement  between  -eKS,i  and  li  is  found  for  the  B3P86  calculations,  the 
magnitude  of  the  average  difference  (-eKS.i  -  Ii)ave  generally  being  less  than  3.0  ev. 
This  is  observed  for  both  the  B3P86/6-31+G**  and  the  B3P86/6-31 1-H-G(3df,3pd) 
calculations;  thus  the  considerable  increase  in  the  basis  set  has  little  effect. 

(3)  Very  similar  results  are  obtained  with  both  the  B3PW91  and  the  B3LYP 
exchange/correlation  combinations;  their  average  differences,  (-eRS.i  -  Ii)ave.  are 
consistently  about  0.6  ev  greater  in  magnitude  than  the  B3P86  for  the  same  molecule. 
Finally,  the  BLYP  results  are  worse  by  approximately  an  additional  1.7  ev. 


presumably  because  the  hybrid  (B3)  procedures  explicitly  include  a  Hartree-Fock 
exchange  contribution  [15, 17]. 

In  Figure  1  is  shown  the  remarkably  parallel  behavior  of  the  average  differences 
obtained  with  the  various  functional  combinations.  This  is  particularly  notable  in  view  of 
the  uncertainties  associated  with  some  of  the  experimental  ionization  potentials.  Figure  1 
suggests  that  the  several  computed  magnitudes  of  (-Eks,!  “  Ii)ave  for  each  molecule  may 

be  related  to  some  property  of  that  molecule.  Recent  work  has  drawn  attention  to  the 
quantity  0.5  (Imjjj  -  A)  in  relation  to  the  asymptotic  behavior  of  the  Kohn-Sham  effective 

potential  [13, 14],  where  A  is  the  molecule’s  electron  affinity.  However  we  have  thus  far 
found  no  similarity  between  the  variations  of  0.5  (Imin  -  A)  and  the  several 
~  Ii)ave* 

Finally,  in  view  of  the  interest  in  the  quantity  0.5  (Iju|jj  -  A)  as  well  as  in  the 
difference  between  the  Kohn-Sham  highest-occupied  and  lowest-unoccupied  orbital 
energies,  eKS,hoo  and  eRS.luo  [9, 1 1],  Figure  2  shows  the  variation  of  0.5(eKS,Iuo  - 
£KS,hoo)  among  the  different  molecules,  for  each  exchange/correlation  combination.  The 
most  striking  feature  is  the  near-uniformity  of  the  results,  other  than  the  BLYP. 

Conclusions 

The  important  points  coming  out  of  this  work  are  the  following: 

(1)  The  magnitudes  of  the  Kohn-Sham  orbital  energies  differ  significantly  from 
experimental  ionization  potentials.  However  the  amount  by  which  they  differ,  for  a 
given  molecule  and  exchange/correlation  combination,  is  usually  roughly  the  same 
for  all  of  the  valence  orbitals. 

(2)  With  the  present  exchange/correlation  functionals,  the  energy  of  the  highest 
occupied  Kohn-Sham  orbital  does  not  provide  a  better  estimate  of  the 
corresponding  ionization  potential  than  do  the  other  valence  orbital  energies. 


(3)  The  average  differences  between  the  magnitudes  of  the  Kohn-Sham  orbital  energies 
and  the  experimental  ionization  potentials  vary  from  molecule  to  molecule. 

However  they  are  always  least  for  the  B3P86  calculations,  about  0.6  ev  greater  for 
the  B3PW91  and  B3LYP,  and  approximately  an  additional  1.7  ev  larger  for  the 
BLYP. 

(4)  These  remarkable  consistencies  suggest  that  the  differences  may  be  related  to  some 
property  of  the  molecule. 
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Table  1.  Differences  between  computed  orbital  energies  and  experimental  ionization  potentials,  in  ev.^ _ 

Molecule  Molecular  -eHF,i  -  li  -CKS.i  -  li 

orbital,  (pi  STO-3G  3-21G  6-3 IG*  6-3 1+G**  B3PW911’  BLYP*’  B3LYPl’  B3P86*’  B3P86c 
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(continued) 


Table  1.  Differences  between  computed  orbital  energies  and  experimental  ionization  potentials,  in  ev  (continued),^ _ 

Molecule  Molecular  -eHF,i  -  li  -eKS,i  -  li 

_ orbital,  cpj  STO-3G  3-2 IG  6-31G*  6-31+G**  B3PW9lb  BLYPb  B3LYPb  B3P86b  B3P86c 
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